Catch-Disperse-Release Readout for Superconducting Qubits 
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We analyze single-shot readout for superconducting qubits via controlled catch, dispersion, and release of 
a microwave field. A tunable coupler is used to decouple the microwave resonator from the transmission line 
during the dispersive qubit-resonator interaction, thus circumventing damping from the Purcell effect. We show 
that if the qubit frequency tuning is sufficiently adiabatic, a fast high-fidelity qubit readout is possible even in the 
strongly nonlinear dispersive regime. Interestingly, the Jaynes-Cummings nonlinearity leads to the quadrature 
squeezing of the resonator field below the standard quantum limit, resulting in a significant decrease of the 
measurement enor. 

PACS numbers: 03.67.Lx, 03.65.Yz, 42.50.Pq. 85.25.Cp 



Introduction.-Fast high-fidelity qubit readout plays an im- 
portant role in quantum information processing. For super- 
conducting qubits various nonline^processes have been used 
to realize a single-shot readout |1]J-|6|]. Linear dispersive read- 
out in the circuit quantum electrodynamics (cQED) setup 
10, S became sufficiently sensitive for the single-shot qubit 
measurement only recently JsKl^l, with development of near- 
quantum-limited superconducting parametric amplifiers |j9|- 
1 in . In particular, readout fidelity of 94% for flux qubits (|91] 
and 97% for transmon qubits ||lo!l has been realized (see also 
1 1211). With increasing coherence time of superconducting 
qubits into 10-100 ^s range iflsllTill . fast high-fidelity readout 
becomes practically important, for example, for reaching the 
threshold of quantum error correction codes lllSll . for which 
the desired readout time is less than 100 ns, with fidelity above 
99%. 

A significant source of error in the currently available 
cQED readout schemes is the Purcell effect lll6ll — the cavity- 
induced relaxation of the qubit due to the always-on coupling 
between the resonator and the outgoing transmission line. The 
Purcell effect can be reduced by increasing the qubit-resonator 
detuning; however, this reduces the dispersive interaction and 
increases measurement time. Several proposals to overcome 
the Purcell effect have been put forward, including the use of 
the Purcell filter lll7ll and the use of a Purcell-protected qubit 
iflill . Here we propose and analyze a cQED scheme which 
avoids the Purcell effect altogether by decoupling the res- 
onator from the transmission line during the dispersive qubit- 
resonator interaction. 

Main idea and reiM/fi. -Similar to the standard cQED mea- 
surement |7-ici, in our method (Fig. 1) the qubit state af- 
fects the dispersive shift of the resonator frequency, that in 
turn changes the phase of the microwave field in the resonator, 
which is then measured via homodyne detection. However, 
instead of measuring continuously, we perform a sequence of 
three operations: "catch", "disperse", and "release" of the mi- 
crowave field. During the first two stages a tunable coupler 
ifigl. Eol] decouples the outgoing transmission line from the 
resonator This automatically eliminates the problems associ- 
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FIG. I: (color online), a) Schematic of the measurement setup. The 
radio frequency (RF) source produces a microwave pulse, which pop- 
ulates the resonator via a small capacitor Cm- The resonator photons 
then interacts with a capacitively (Cg) coupled qubit. The interac- 
tion with the outgoing transmission line is controlled by a tunable 
coupler, which releases photons at the end of the procedure. The re- 
leased field is then amplified and mixed with the local oscillator (LO) 
signal to be measured via homodyne detection, b) The RF pulse B{t) 
(blue curve) and varying qubit frequency cjq(t) (red curve), with ap- 
proximate indication of the "catch", "disperse", and "release" stages. 
Dashed lines show the resonator frequency uii and initial/final qubit 



frequency uq; A. — 



is the detuning at the "disperse" stage. 



ated with the Purcell effect, as coupling to the incoming mi- 
crowave line can be made very small COll . 

During the "catch" stage, the initially empty resonator is 
driven by a microwave pulse and populated with ~^10 pho- 
tons. At this stage the qubit is far detuned from the resonator 
[Fig. 1(b)], which makes the dispersive coupling negligible 
and allows the creation of an almost-perfect coherent state in 
the resonator At the next "disperse" stage of the measure- 
ment, the qubit frequency is adiabatically tuned closer to the 
resonator frequency to produce a strong qubit-resonator in- 
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teraction (it may even be pushed into the nonlinear regime). 
During this interaction, the resonator field amplitudes (Acff) 
associated with the initial qubit states |0) and |1) rapidly ac- 
cumulate additional phases and separate in the complex phase 
plane [see Fig. |2ja)]. Finally, at the last "release" stage of 
the measurement, after the qubit frequency is again detuned 
from the resonator, the resonator photons are released into the 
outgoing transmission line. The signal is subsequently ampli- 
fied (by a phase-sensitive parametric amplifier) and sent to the 
mixer where the homodyne detection is performed. 

With realistic parameters for superconducting qubit tech- 
nology, we numerically show that the measurement of 30-40 
ns duration can be realized with an error below lO^'^, neglect- 
ing the intrinsic qubit decoherence. The latter assumption re- 
quires the qubit coherence time to be over 40 ^s, which is 
already possible experimentally lll4ll . It is interesting that be- 
cause of the interaction nonlinearity 11211 l22ll . increasing the 
microwave field beyond ~ 10 photons only slightly reduces the 
measurement time. The nonlinearity also gives rise to about 
~50% squeezing of the microwave field, which provides an 
order-of-magnitude reduction of the measurement error 

The model.-We consider a superconducting phase or trans- 
mon qubit capacitively coupled to a microwave resonator [Fig. 
[Ha)]. For simplicity we start with considering a two-level 
qubit (the third level will be included later) and describe the 
system by the Jaynes-Cummings (JC) Hamiltonian 0] with a 
microwave drive (we use h ^ 1) 



H = cjq(t)CT+CT_ + Lo^a^a + g{aa+ + a-o)) 



(1) 



where Wq(i) and Wj. are, respectively, the qubit and the res- 
onator frequencies, a± are the rasing and lowering operators 
for the qubit, a (a^) is the annihilation (creation) operator for 
the resonator photons, g (assumed real) is the qubit-resonator 
coupling, B{t) and uj are the effective amplitude and the fre- 
quency of the microwave drive, respectively. In this work we 
assume uj = uj-c- 

For the microwave drive B{t) and the qubit frequency cjg {t) 
[Fig. 1(b)] we use Gaussian-smoothed step-functions: B{t) = 
0.5Bo{Erf[(t - tB)/\/2<TB] - Erf[(t - - m)|^/2a^]} and 
cjq(i) = Wo + 0.5(Ao - A){Erf[(t - t^)lV2a^] - Erf[(t - 
tqo)/-\/2(Tqe]}, whcre Ib, t-Q + tb, tq, and iqc are the centers 
of the front/end ramps, and ctb, cTq, and cTqo are the corre- 
sponding standard deviations. In numerical simulations we 
use (Tb = CTqc = 1 ns (typical experimental value for a short 
ramp) while we use longer cTq to make the qubit front ramp 
more adiabatic. Other fixed parameters are: g/2'K = 30 MHz, 
Tb = 1 ns, t-Q ^ 'i ns, uj^jl-K = 7 GHz, and uj^II-k = 6 GHz, 
so that initial/final detuning Aq = w, — wq is 1 GHz, while the 
disperse-stage detuning A is varied. The measurement starts 
at t = and ends at tf = tqc + la^c when the field is quickly 
released 123 1. 

Simplified analysis-Let us first consider a simple dis- 
persive scenario at large qubit-resonator detuning, |A| ^ 
g\^n + 1, where n is the average number of photons in 
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FIG. 2: (color online), (a) Evolution in time of the effective field am- 
plitude Acff on the complex phase plane for initial qubit states |0) and 
|1), computed numerically. The horizontal line emphasizes slight 
asymmetry. The dots indicate time moments t — 0, 5, 10, 15, 20, 
25, and 30 ns. (b) Corresponding probability distributions Po{x^p) 
and Pi(x^) for measurement (at t — t{) of the optimum quadrature 
x^p. Side bumps of Po and Pi are due to non-adiabaticity. We used 



g/2n = 30 MHz, A/27r = 50 MHz, |A,„r = 
MHz, Tb = 1 ns), — 3 ns, tq — 3.25 ns, tq 
and ti = 32 ns. 



9 {Bo/2n = 497.4 
= 30 ns, tB ~ S ns, 



the resonator In this case, the system is described by the 
usual dispersive Hamiltonian 101 Hd = {oJo + g^/A)(T^/2 + 
{ujr + <Tzg'^ /A)a^a, where az is the Pauli matrix. Af- 
ter the short "catch" stage the system is in a product state 
(a|0) + /3|l))|Ain), where a and /3 are the initial qubit state 
amplitudes and Ain is the amplitude of the coherent resonator 
field, Ain = J B{t)dt (so n = |AinP). Then during the 
"disperse" stage the qubit-resonator state becomes entangled, 
a|0)|Ao(i))+/3|l)|Ai(i)), withAo = Ai„e-*^,Ai = Ai„e*^ 

and0(i)=/„*feVA(t')]rft'- 

The distinguishablity of the two resonator states depends on 

their separation \6X\ = |Ai — Ao| — 2|Ain| sin (see numer- 
ical results in Fig. m. The released coherent states are mea- 
sured via the homodyne detection using the optimal quadra- 
ture connecting Aq and Ai, i.e. corresponding to the angle 
ip = arg(Ai — Ao). We rescale the measurement results to 
the dimensionless field quadrature = (ae~"^ + a^^e*'^)/2, 
which corresponds to the iy9-angle axis in the phase space of 
Fig. |2ja). In resolving the two coherent states, we are essen- 
tially distinguishing two Gaussian probability distributions, 
Pa{x^) and Pi{xip), centered at il^Alo-coh with dcoh = 1/2 
being the coherent-state width (standard deviation) for both 
distributions. Then the measurement error has a simple form 



E 



min(Po, -Pi) dx^ 



1 -Erf(|<5A|, 



(2) 



where rj ~ ?/coif;amp is the detection efficiency |24j| . which 
includes the collection efficiency t^coI and quantum efficiency 
of the amplifier r^amp- Unless mentioned otherwise, we as- 
sume rj = 1, which corresponds to a quantum-limited phase- 
sensitive amplifier (for a phase-preserving amplifier 77 < 1/2). 

Full analysis.-ln general the JC qubit-resonator interac- 
tion (dl) is non- linear for |Ai„|2 > A^/4g^ H and the res- 
onator states are not coherent. The measurement error E is 
still given by the first part of Eq. (|2]), while the probabil- 
ity distributions Pq.i{xip) of the measurement result for the 
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qubit starting in either state |0) or |1) can be calculated in 
the following way. Assuming an instantaneous release of the 
field, we are essentially measuring the operator i^. There- 
fore the probability P{x^) for the ideal detection (77 = 1) 
can be calculated by converting the Fock-space density ma- 
trix pnm describing the resonator field, into the cc^-basis, thus 



obtaining P{x^) = T,nm^nix^)pnrnit)lp*aix^)e- 

where ipn{x) is the standard nth-level wave function of a har- 
monic oscillator For a non-instantaneous release of the mi- 
crowave field the calculation of P{x^) is non-trivial; how- 
ever, since the qubit is already essentially decoupled from 
the resonator, the above result for P{x^p) remains the same 
lEsll for optimal time-weighting of the signal. In the case of a 
non-ideal detection (77 < 1) we should take a convolution of 
the ideal P{x^) with the Gaussian of width \/rj^^ — 1 (Tcoh- 
Calculation of the optimum phase angle ip minimizing the er- 
ror is non-trivial in the general case. For simplicity we still 
use the natural choice ip = arg(Aeff.i — A^ff^o), where the 
effective amplitude of the resonator field ll26ll is defined by 
Aeff = J2n V^Pn,n-i- The field density matrix pnm is calcu- 
lated numerically using the Hamiltonian ([T]i and then tracing 
over the qubit. 

Extensive numerical simulations allowed us to identify two 
main contributions to the measurement error E in our scheme. 
The first contribution is due to the insufficient separation of 
the final resonator states |Aoff.i) and |Aoff,o), as described 
above. However, there are two important differences from 
the simplified analysis: the JC nonlinearity may dramatically 
change \SX\ and it also produces a self-developing squeez- 
ing of the resonator states in the quadrature .x^, significantly 
decreasing the error compared with Eq. dH (both effects are 
discussed in more detail later). The second contribution to 
the measurement error is due to the nonadiabaticity of the 
front ramp of the qubit pulse ujq{t), which leads to the pop- 
ulation of "wrong" levels in the eigenbasis. This gives rise 
to the side peaks ("bumps") in the probability distributions 
Po,i{xip), as can be seen in Fig.Ulb) (notice their similarity to 
the experimental results IsKI^I, though the mechanism is dif- 
ferent). During the dispersion stage these bumps move in the 
"wrong" direction, halting the exponential decrease in the er- 
ror, and thus causing the error to saturate. The nonadiabaticity 
at the rear ramp of ujq{t) is not important because the moving 
bumps do not have enough time to develop. Therefore the rear 
ramp can be steep, while the front ramp should be sufficiently 
smooth [Fig. 1(a)] to minimize the error 

Now let us discuss the effect of nonlinearity (when | Ainp > 
A^/4^^) on the evolution of Acff,o and Aoff,i during the 
disperse stage. Since the RF drive is turned off, the in- 
teraction described by the Hamiltonian ^ occurs only be- 
tween the pairs of states |0,n) and \l,n — 1) of the JC 
ladder Therefore, if the front ramp of the qubit pulse 
is adiabatic, the pairs of the JC eigenstates evolve only 
by accumulating their respective phases while maintain- 
ing their populations. Then for the qubit initial state |0), 
the qubit-resonator wavefunction evolves approximately as 
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FIG. 3: (color online). Optimized measurement error E vs measure- 
ment time ti (optimization is performed over A, Uq, and tq). The mi- 
crowave pulses correspond to mean photon number | Ainp — 6, 9, 12, 
and 15. The lines are guides for the eye. 



the overbar denotes the (dressed) eigenstate and 0o.n(i) = 
/j*^ (ii'[^A(t')^ + A(t')]/2 is the accumulated phase, 

with t-o = ie + t'b/S being the center of the i?(<)-pulse, 
which is crudely the start of the dispersion. Similarly, if 
the qubit starts in state |1) (following the ideology of Ref. 
1 2711 . we then use |10) as the initial state), the state evolves 
as iVi(i)) - e-l^i»l'/2X;„(A[;,/%/^)e''^i'"(*)|T;7I), where 
^lAt) = !l dt'WA{tr + ^g^n + 1) - A{t')]/2. Us- 
ing the above definition of Acff and assuming |AinP ^ 1 we 
derive an approximate formula 



Aeff.o = Ain exp 



■.dt' 



(3) 



The corresponding expression for Aeff,i can be obtained by re- 
placing —i with i and |AinP with |AinP + 1- These formulas 
agree well with our numerical results. In particular, they ex- 
plain why Aeff rotates slightly faster than Aeff 1, as seen in 
Fig. Ha). 

Equation Q shows that a decrease in detuning leads to an 
increase in the rotation speed of Aeff. However, in the strongly 
nonlinear regime |Ainp S> the angular speed sat- 

urates at (i(arg(Aeff,o/i))/cJt = T§/2|Ain|. Thus, the rate at 
which the Aeff,i and Agff.o separate is limited by 



d\SX\/dt < 



(4) 



\Mt)) 



-|A,„|V2 



-100, 



O.n), where 



which does not depend on |Ain|. This means that the mea- 
surement time should not improve much with increasing the 
mean number of photons |Ainp in the resonator, as long as it 
is sufficient for distinguishing the states with a desired fidelity 
(crudely, |Ainp > 7/r] for E < IQ-*). 

Results of numerical optimization.-Figme^sho'ws the re- 
sults of a three-parameter optimization of the measurement 
error E for several values of the average number of photons 
in the resonator, |Ainp (assuming jy = 1). The optimization 
parameters are the qubit-resonator detuning A, the width dq, 
and the center tq of the qubit front ramp. We see that for 9 
photons in the resonator the error of 10~^ can be achieved 
with 30 ns measurement duration, excluding time to release 
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;; = 0.5, withlevel|2> 
7/ = 0.5, without level 1 2) 
;;= 1.0, with level 1 2> 
r]= 1.0, without level 1 2) 
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FIG. 4: (color online). Optimized error vs measurement time t{ for 
|Ai„|^ — 9 and quantum efficiencies rj = 1 or rj — 1/2 (e.g. for a 
phase-preserving amplifier), taking into account the qubit level |2) 
(with anharmonicity ^/27r = 200 MHz) or assuming a two-level 
qubit. 



and measure the field. The optimum parameters in this case 
are: A/27r = 60 MHz, cr, 4.20 ns, and = 3.25 ns [this 
is a strongly nonlinear regime: | Ainp/(A^/4§^) = 9]. As 
expected from the above discussion, increasing the mean pho- 
ton number to 12 and 15 shortens the measurement time only 
slightly (by 1 ns and 2 ns, keeping the same error). The dotted 
blue curve in Fig. |4] shows the optimized error for |Ainp = 9 
and imperfect quantum efficiency rj = 1/2. As we see, the 
measurement time for the error level of 10"'* increases to 40 
ns, while the error of 10^^ is achieved at if = 32 ns. 

So far, we considered the two-level model for the qubit. 
However, real superconducting qubits are only slightly anhar- 
monic oscillators, so the effect of the next excited level |2) 
is often important. It is straightforward to include the level 
1 2) into the Hamiltonian ([T]i by replacing its first term with 
a;q|l)(l| + (2tjq — ^)|2)(2|, where A is the anharmonicity. 
The dispersion can then be understood as due to repulsion 
of three eigenstates: |0,n), |l,n— 1), and |2,n — 2). As 
the result, Acff ,o rotates on the phase plane faster than in the 
two-level approximation, while Acff.i rotates slower (some- 
times even in the opposite direction). The Supplemental Ma- 
terial ll28ll illustrates evolution of the resonator Wigner func- 
tion corresponding to initial qubit states |0) and |1). The de- 
tailed description of the effect of level |2) is beyond the scope 
of this paper Here we only show the optimized error for 
A/27r ~ 200 MHz (a typical value for transmon and phase 
qubits), |AinP = 9 and 77 = 1 or 77 = 1/2 as the red lines in 
Fig- El For the error of lO^'^ the measurement time tf is 31 ns 
and 39 ns, respectively. 

The presented results are for the coupling g/2n = 30 MHz. 
While the dispersion rate scales as g in the strongly nonlinear 
regime [see Eq. (01)], the scaling of the overall duration of the 
measurement process is slower than because of signifi- 
cant time needed at the initial and final stages [see Fig. 1(b)]. 

In this work we are not interested in how close our read- 
out is to the quantum non-demolition (QND) measurement 
I 29n . Notice that in the proposed implementation of the sur- 
face code ifisll the measured qubits are reset, so the QND- 
ness is not expected to be important. For the results presented 
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FIG. 5: (color online), (a) Time-evolution of the quadrature squeez- 
ing (the qubit is initially in state |0)). (b) Measurement error vs |(5A 
obtained numerically (solid lines) and using Eq. (|2j (dashed line); 
the evolution stops at 98 ns. Here g/27r = 30 MHz, |A,„p = 9, 
3.25 ns. 



in Fig. O the non-QNDness (probability that the initial states 
1 00) and |10) are changed after the procedure, including the 
release stage) is crudely about 5%, which is mainly due to 
non-adiabaticity of the rear ramp of the qubit pulse. It is pos- 
sible to decrease the non-QNDness significantly by using a 
larger rear ramp width <Tqo and correspondingly increasing 
the overall duration by several nanoseconds. However, the 
non-QNDness cannot be reduced below the level of few times 
(g/Ao)^, essentially because of the PurceU effect during the 
release stage. 

Squeezing.-The JC nonlinearity causes the quadrature 
squeezing of the microwave field. To quantify the squeez- 
ing, we calculate the variance Ax'^ = (x^) — {x^)"^, which 
can be put in the form Aa;^ = 1/4 + {a^a)/2 - |(a)p/2 + 
Rc[((a2) - (a)2)e-2iv]/2. For a coherent field Axl = 1/4, 



4Az2 



thus the state is squeezed Il26ll when 4Aa;^ < 1. Notice that 
the degree of squeezing depends on the choice of ip. However, 
in order to see the effect of squeezing on measurement error, 
we compute the squeezing along the measurement direction, 
V? = arg(Aeff,i - Aeff,o). 

Figure [3 a) shows evolution of the squeezing parameter 
when the initial qubit state is |0) (a similar result is 
obtained for the initial state |1); we again use the two-level 
approximation and i] = 1). Notice that at first the field 
stays coherent, which is due to the linearity of the qubit- 
resonator interaction at large detuning. Later on, however, the 
interaction becomes nonlinear due to decreased detuning and 
leads to quadrature squeezing reaching the level of ~^50% for 
A/27r < 100 MHz (see [28)1 for the Wigner function evo- 
lution). Figure |5lb) shows the measurement error as a func- 
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tion of |(5A| calculated numerically (solid curves) and analyti- 
cally, based on Eq. ^ (dashed curve). As expected, with the 
squeezing developing, the error becomes significantly smaller 
than the analytical prediction, for instance, up to a factor of 
30 for A/27r = 250 MHz. Notice also that the error shown in 
Fig-Etb) saturates in spite of increasing separation \6X\. This 
is because of the non-adiabatic error discussed above. 

Conclusion.-We analyzed a fast high-fidelity readout for 
superconducting qubits in a cQED architecture using the con- 
trolled catch, dispersion, and release of the microwave pho- 
tons. This readout uses a tunable coupler to decouple the res- 
onator from the transmission line during the dispersion stage 
of the measurement, thus avoiding the Purcell effect. Our ap- 
proach may also be used as a new tool to beat the standard 
quantum limit via self-developing field squeezing, directly 
measurable using the state-of-the-art parametric amplifiers. 
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